Atom-wave diffraction between the Raman-Nath and the Bragg regime: 
Effective Rabi frequency, losses, and phase shifts. 
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We present an analytic theory of the diffraction of (matter) waves by a lattice in the "quasi- 
Bragg" regime, by which we mean the transition region between the long-interaction Bragg and 
"channelling" regimes and the short-interaction Raman-Nath regime. The Schrodinger equation is 
solved by adiabatic expansion, using the conventional adiabatic approximation as a starting point, 
and re-inserting the result into the Schrodinger equation to yield a second order correction. Closed 
expressions for arbitrary pulse shapes and diffraction orders are obtained and the losses of the 
population to output states otherwise forbidden by the Bragg condition are derived. We consider 
the phase shift due to couplings of the desired output to these states that depends on the interaction 
strength and duration and show how these can be kept negligible by a choice of smooth (e.g., 
Gaussian) envelope functions even in situations that substantially violate the adiabaticity condition. 
We also give an efficient method for calculating the effective Rabi frequency (which is related to the 
eigenvalues of Mathieu functions) in the quasi-Bragg regime. 

PACS numbers: 03.75.Be; 32.80.Lg; 32.80.Wr; 03.75.Dg 



I. INTRODUCTION 

A. Background 

Diffraction by a point scatters light or matter waves 
into all directions. A two-dimensional grating produces 
a few diffraction orders at those angles where the scatter 
from all of the grating adds coherently. Bragg diffrac- 
tion by an infinite three-dimensional lattice can produce 
a single diffraction order, which happens when the scat- 
ter from all layers adds constructively, as described by 
the Bragg condition. When this happens for a higher 
scattering order ("high-order Bragg diffraction") virtu- 
ally all incident radiation can be scattered into this high 
order, in contrast to the two-dimensional case. By quasi- 
Bragg diffraction, we refer to the intermediate regime 
where the infinite lattice assumption is no longer valid 
but approximately true. Using the nomenclature of, 
e.g., [J, this regime is the transition between the short- 
interaction Raman-Nath regime and the long-interaction 
Bragg (weak potential) and "channelling" (strong poten- 
tial) regimes. In this region, the Bragg condition soft- 
ens and there may be significant scattering into other 
than the desired orders. Moreover, couplings between 
the nonzero diffraction orders may lead to phase shifts 
of the diffracted waves [H , which is undesirable in many 
applications. 

In this work, we present an analytic treatment of such 
quasi Bragg scattering. We will find that by prudent 
choice of the scattering potential and its envelope func- 
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tion, behavior very similar to Bragg scattering, in par- 
ticular very low losses and phase shifts, can be obtained 
for scatterers that substantially violate the assumptions 
of the simplified theory. 

Bragg diffraction famously provides us with the basic 
knowledge of the structure of crystals, including proteins. 
It is also important for many technical applications, like 
acousto-optic modulators ( AOMs) [3[ , distributed Bragg 
reflectors (DBR) in diode [4 and fiber lasers as well as 
photonic bandgap crystals @. Moreover, Bragg diffrac- 
tion is a basic method for making surface acoustic wave 
(SAW) filters in radio frequency technology. In atomic 
physics, Bragg diffraction is a special case of the Kapitza- 
Dirac effect fj|, 0, i, i, El ■ 

Bragg scattering is used as a tool for experiments with 
Bose-Einstein condensates (BECs) [ill El Ell • For ex- 
ample, Kozuma et al. [l4j have shown experimentally 
that thirteen subsequent first-order Bragg diffractions of 
a BEC can still have good efficiency. More exotic ap- 
plications include the generation of a collective frictional 
force in an ensemble of atoms enclosed in a cavity, due 
to Bragg scattering of a pump light an a self-organized 
atomic density grating [15], much in the same way as 
stimulated Brillouin scattering by self-organized acoustic 
waves in optical fibers Q. 

Moreover, Bragg diffraction can act as a beam split- 
ter for matter waves 0, 0, EH, El El, El El El 113] • 
The highest order diffraction so far achieved with matter 
waves seems to be by Koolen et al. [2l| , who obtained up 
to eighth-order Bragg diffraction. Atom interferometers 
based on Bragg diffraction include the one by Giltncr 
et al. [H l23j ] who built a Mach-Zehnder atom inter- 
ferometer using up to third order diffraction. Miller et 
al. [Hj] achieved high contrast in a two-pulse geometry 
with first-order diffraction and a sufficiently short time 
between pulses. Torii et al. [25| have used first order 
Bragg diffraction in a Mach-Zehnder geometry with a 
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Bose-Einstein condensate. In addition, Rasel et al. [261 ] 
have built a Mach-Zchnder atomic-beam interferometer 
based on Raman-Nath scattering. 

More generally, atom interferometers can be used for 
measurements of atomic properties [13, [H, HU, the lo- 
cal gravitational acceleration [301] . t he gravity gradient, 
Newton's gravitational constant [3l|, tests of the equiv- 
alence principle [32| and the fine-structure constant via 
jt/m I33l . l34l .l35j]. For planned experiments in space, see 

While not all the atom interferometers just cited use 
Bragg diffraction, high-order Bragg diffraction offers sev- 
eral interesting possibilities for atom interferometers: (i) 
it makes the atom interact with 2n photons at once, 
which may increase the sensitivity of the interferometer 
by a factor of n 2 [35| relative to the 2— photon tran- 
sitions. (We note here that other possibilities exist for 
using high-order transitions in atom interferometers, like 
applying multiple low-order pulses [38|, operation in the 
Raman-Nath regime [3{| or the magneto-optical beam 
splitter [13]) (ii) Since Bragg diffraction theoretically al- 
lows coherent momentum transfer with an efficiency close 
to one, it allows the insertion of many 7r-pulses for addi- 
tional momentum transfer, which increases the signal in 
photon recoil measurements 35] (up to N = 30 7r-pulses 
based on two-photon adiabatic transfer were used in [33j , 
transferring 60 photon momenta; if these had been 5-th 
order Bragg pulses, they would have transferred 300). 
(iii) If losses can be neglected, Bragg diffraction is ba- 
sically a transition in a 2-level system. Thus, many 
of the techniques developed for standard beam splitters 
based on Raman transitions can be taken over. For ex- 
ample, several beam splitters addressing different veloc- 
^roups respectively can be performed simultaneously 
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Overview of the existing theory 



A summary of the material that is the basis of this 
work can be found in the textbook by Meystre [l6l ]. 

A lot of attention has been paid on the theory of the 
long-interaction time (channelling or Bragg) regimes on 
the one hand and the Raman-Nath regime on the other 
hand. Keller et al. [l[ give a brief account of the most 
important results. They have been derived using vari- 
ous formalisms: Bcrman and Bian |9( use a pump-probe 
spectroscopy picture, focussing on applications as beam 
splitters in atom interferometers. The phase-shift of the 
diffraction process has been studied by Biichner et al. Q , 
limited to first and second order diffraction. Giltner et 
al. [13, HH] have reported an atom interferometer based 
on Bragg diffraction of up to third order and give the 
effective Rabi frequency in the long-interaction regime, 
our E g. (1151) . A similar derivation was given by Gupta 
et al. 

While most of this work (as well as ours) is con- 
cerned with on-resonance transitions, Diirr and Rempe 



[42| have considered the acceptance angle (i.e., linewidth) 
of diffraction. They restrict attention to the case of 
square-envelope pulses. Wu et al. [H[ give a numerical 
study of the latter, again for the case of (a pair of) square 
pulses. Stenger et al. [43| describe the line shape of 
Bragg diffraction with Bose-Einstein condensates within 
the Bragg regime. 

More general types of diffraction have been studied, 
such as using chirped laser frequencies [13] in the adia- 
batic regime where the chirp is sufficiently slow. This 
regime is similar to Bloch oscillations [44| . Band et 
al. [I]! have considered the loss due to atom-atom in- 
teractions, which is relevant in experiments with Bose- 
Einstein condensates (BECs). Blackie and Ballagh [46[ 
explore the use of Bragg diffraction in probing vortices 
in BECs. 

On the other hand, the picture is still incomplete in 
the quasi-Bragg regime we are concerned with. Diirr and 
Rempe have analytically calculated corrections to the ef- 
fective Rabi frequency for large potential depth [47[ in 
the case of second-order scattering with a square enve- 
lope function. Champenois et al. [481 ] consider both the 
Raman-Nath and the Bragg regime using a Bloch-state 
approach; for the Bragg regime, however, their treatment 
is restricted to the first diffraction order and square en- 
velopes. Borde and Lammerzahl [l9| give and exhaustive 
treatment of square-envelope scattering that is based on 
matching the boundary conditions at the beginning and 
end of the pulse. This method, unfortunately, cannot 
readily be generalized to smooth envelope functions. 

The Mathieu equation formalism is a powerful tool 
for studying diffraction with arbitrary potential depth 
and interaction times, as demonstrated by Home, Jex, 
and Zeilinger [49|. However, like the work of Borde and 
Lammerzahl, the Mathieu equation formalism assumes 
constant envelopes and cannot readily be generalized to 
smooth envelope functions. As we shall see (and has 
already been pointed out, see, e.g., Q), smooth enve- 
lope functions are of particular interest because they al- 
low high-efficiency scattering into a single order even in 
the quasi-Bragg regime. The mathematical properties of 
the Mathieu functions themselves have been explored by 
many workers. Of relevance for this work is the power- 
series expansion of Mathieu functions reported by Kokko- 
rakis and Roumelotis [501 ■ 



C. Motivation 

In the experimental applications, it is often desirable 
to make the interaction time as short as compatible with 
certain requirements on the efficiency and parasitic phase 
shifts. This can be achieved by operating in the quasi- 
Bragg regime. For example, in atomic physics, long in- 
teraction times increase losses due to single-photon ex- 
citation and also systematic effects in atom interferome- 
ters. For fifth order Bragg scattering of cesium atoms in 
a standing light wave, satisfying the adiabaticity crite- 
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rion requires interaction times 3> 0.4s, see Sec. HE 21 
This exceeds the time available in experiments under 
free-fall conditions, and would give rise to huge losses 
by single-photon excitation. Moreover, the very sharp 
Bragg condition in the pure Bragg regime means that 
scattering happens only if the incident waves are within 
a very narrow of the velocity distribution of a thermal 
sample, which may mean that a large fraction will not be 
scattered at all. However, operation in the quasi-Bragg 
regime requires a theoretical calculation of the losses 
and phase-shifts encountered, and strategies to minimize 
them. This is best done by analytic equations for these 
parameters, which allow one to easily see which parame- 
ters have to take which values. 

Unfortunately in this regime the population of the 
diffraction orders is particularly difficult to calculate [l[ 
and so the existing analytic theory of the quasi-Bragg 
regime is restricted to square envelope functions or low 
orders, or both. Moreover, even then, using the known 
formalisms it is hard to obtain power-series approxima- 
tions of parameters like the effective Rabi frequency in 
terms of the interaction strength. 

The aim of this paper is to present an analytic the- 
ory of the quasi-Bragg regime that allows the treatment 
of arbitrary scattering orders and envelope functions of 
the scattering potential. To do so, we develop a system- 
atic way of obtaining more and more accurate solutions 
of the Schrodinger equation that starts from the usual 
adiabatic approximation. This allows us to calculate the 
population of the diffraction orders, including losses to 
unwanted outputs and phase shifts. We can thus specify 
the minimum interaction time and the maximum inter- 
action strength that yield losses which are below a given 
level, and consider the influence of the pulse shape. It 
will turn out that efficient scattering can be maintained 
with interaction times that substantially violate the adi- 
abaticity criterion, in agreement with experiments 
For example, we show that fifth order scattering of Cs 
atoms still has negligible losses and phase-shifts for in- 
teraction times on the order of 10 /js if the pulse shape 
of the light is appropriately chosen. This minimum in- 
teraction time even decreases for higher scattering order. 
We will restrict attention to on-resonant Bragg diffrac- 
tion, neglecting any initial velocity spread that the atoms 
may have. This can certainly be a good assumption for 
experiments using Bose-Einstein condensed atoms, but 
also for a much wider class of experiments: The min- 
imum interaction time that will still lead to low losses 
will turn out to be roughly given by l/(nu> r ), where to r 
is the recoil frequency and n the Bragg diffraction order, 
see Eq. (|74l) . Such short pulses have a Fourier linewidth 
that is on the order of noj r . This means that a velocity 
spread on the order of the recoil velocity cannot be re- 
solved, even though the velocity selectivity increases ~ n 
due to the multiple photons scattered. The preparation 
of atoms having a velocity spread of 1/100-1/10 of the 
recoil velocity is standard practice in atomic fountains. 

While our theory will be stated in the language of 



atomic physics, with an eye on applications in atom in- 
tcrferomctry, it can be adapted to Bragg diffraction in 
all fields of physics. 

D. Outline 

This paper is organized as follows: In Sec. II El we de- 
scribe our basic Hamiltonian and the conventional theory 
for the Bragg and the Raman-Nath regime. In Sec. II F[ 
an exact solution for rectangular envelope functions will 
be presented that uses Mathieu functions. In Sec. [Til 
we find a general form for the corrections of the effective 
Rabi frequency which is valid for arbitrary scattering or- 
ders. In Sec. Hill wc present our method for calculating 
the phase shifts and losses, adiabatic expansion. In Sec- 
tions [IV] and [Vj wc consider square and Gaussian enve- 
lope functions and give a practical example of high-order 
Bragg scattering of Cs atoms. 

E. Problem 

In the remainder of this Sec. I, we define the basic 
problem and review the basic theory of the adiabatic and 
the Raman-Nath case (as described, e.g., in and of 
the Mathieu equation approach. This is to define the 
notation and for the reader's convenience. The reader 
already familiar with this may want to proceed to the 
following sections, which describe the new results of this 
paper. 

Consider scattering of an atom of mass M by a stand- 
ing wave of light along the z direction. Ignoring effects 
of spontaneous emission, the Hamiltonian describing the 
interaction of the atoms with the standing wave having 
a wavenumber k is (in a frame rotating at the laser fre- 
quency u)) 

2 

H=^-h5\e)(e\ + hCl (t)cos(kz)(\e)(g\+H.c.) , (1) 

where [z,p] = ih and S is the detuning. The Rabi fre- 
quency fi may in general be time-dependent. For the 
purpose of this introduction, we will assume it to be con- 
stant. In the later sections of this paper, we shall be 
interested in the effects of different pulse shapes, how- 
ever. Substituting 

\4>(t))=e(z,t)\e)+g(z,t)\g) (2) 

into the Schrodinger equation yields the coupled differ- 
ential equations 

P 2 

ihe(z,t) = — ~£ e \ z i v + f&lo cos(k z)g(z,t) — hSe(z,t) , 
P 2 

ihg(z,t) — ——g(z,i)+hSlocos(kz)e(z,t), (3) 

where the dot denotes the time derivative. For d large 
compared to the linewidth of the excited state (and thus 
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also 6 ^> Qo,u r ) and the atoms initially in the ground 
state, we can adiabatically eliminate the excited state: 



ihg(z, t) 



h 2 d 2 g(z,t) 
2M dz 2 



^cos 2 (kz)g(z,t). (4) 



This equation with its periodic potential is invariant un- 
der a translation by an integer multiple of fc _1 . Apply- 
ing the Bloch theorem, we can look for solutions having 
constant quasi-momcntum; in particular, we can restrict 
attention to the case of vanishing quasi-momentum. For 
constant Hq, this is a Mathieu equation for which exact 
solutions are known; this formalism will be described in 
Sec. II Fl and developed further in Sec|TTJ If we let 



]( z ,t)= ^2 9m(t)( 



(5) 



and use cos 2 (kz) = 1/2 + l/4(e 2lkz + e 2tkz ) i we obtain 

oo oo 

ih V g m (t)e mkz =h V [(cj r m 2 +n)g m 



7TL— — OO 



m= — oo 



+ (n/2)(g m+2 +g m ^)]e mkz (6) 
where we have introduced the two-photon Rabi frequency 



O 2 

fi = 

26 



and the recoil frequency 



hk 2 

2M 



(7) 



(8) 



This can only hold if for all m 

ihg m = h(u> r m 2 + Q)g m + H(ft/2)(g m+2 + g m -2) ■ (9) 

[Since this equation couples only odd or even momentum 
states, respectively, we can and will look for solutions 
that have either the odd or even terms zero. In view 
of this the use of both even and odd indices may seem 
unnecessary, but will have advantages when we consider 
Bragg diffraction.] The theoretical description of Bragg 
diffraction is relatively simple in the short-interaction 
limit (the Raman-Nath regime) and in the case of an 
infinite scatterer, the Bragg regime. 

1. Raman-Nath Regime 

The Raman-Nath regime is defined as the case of very 
short interaction time, so that the kinetic energy term is 
negligible against the resulting energy uncertainty. Equa- 
tion © reduces to 



ihg m = h(fl/2)(g m+2 + gm-2) , 



(10) 



which we have simplified by shifting the energy scale by 
— HQ. Since these equations only couple states which dif- 
fer by an even multiple of the momentum hk, we can 



restrict attention to even indices 2m. They can be satis- 
fied by Bessel functions: 



92r. 



H) m J m (Qt). 



(11) 



At t = 0, this solution has all atoms in the zero momen- 
tum state go = l,<?2m^o — 0. For t 7^ 0, the probability 
to find the atom to have a transverse momentum 2mhk 
is P2m(t) = J m {Qt)- The Raman-Nath approximation 
holds provided that t -C 1 / ^/2Qu r (because then a high 
energy uncertainty justifies our neglect of the kinetic en- 
ergy) . Clearly, the transfer efficiency P2m for any particu- 
lar m is limited. For example, the maximum probability 
to find the atom in the ground state after scattering 2 
photons is approximately 0.34. 



2. Bragg Regime 

For the Bragg regime, we take into account the kinetic 
energy term and work in configuration space. We now 
assume initial conditions g_„ = 1 and g m — for m ^ 
— n. To simplify, we subtract a constant offset n 2 hui r +hfl 
from the energy scale. Eq. ([9]) now reads 



ihg-„^ 2 = (4 + An)huj r g- n - 2 + -Ml(g_ n + ff_„_ 4 ) 
ifig-n+2 = (4 - 4n)hw r g- n+ 2 + -ffil(g- n+ 4, + g_ n ) 



ihg^ n+2 k = 4fe(fc - n)fiw r g- n+ 2k 

+ -hfl(g- n+ 2k+2 + g-n+ik-2) 



i%g n -2 = (4- 4n)fiw r g„„ 2 + -hH{g n + 9n-i) 

itl9n = 2^(dn+2 + gn-2) 

ifrg-n+2 = (4 + 4n)huj r g n+ 2 + ^hil(g n+4: + g n ) 



(12) 



Energy conservation will favor transitions from — n — > n, 
if the processes are sufficiently slow. This is the result 
of the adiabatic elimination of the intermediate states 
(fc ^ and k ^ n): If 

|4fc 2 -Ank\hw r > Ml, (13) 

for all < k < n, we can assume that the fcth equation 
is always in equilibrium with g- n +2k ~ 0. Then, for 
example, 



g-n+2 = 



8 {nk — k 2 )hw r ' 



(14) 
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Relations like this can be used to successively eliminate 
all n — 1 intermediate states. With 



a 



9." 



cff 



n—1 

t n „h - 



1 



(Sw,.)"- 1 iJL nk - k 2 {SuJr)*- 1 [(n - l)!] 2 ' 

(15) 



we obtain 



ihg- n = -hn eS g n 



(16) 



(where we have removed a constant light shift term). 
This can be readily solved: 

g- n {t) = cos ^Q e $t , g n (t) = -i sin ^Q e $t . (17) 
For a time- varying Q, 



3 _„(<) = cos ( - / Sl eB (t')dt' 



g n (t) = ~*sin^- J n cS {t')dt'j . (18) 

This is an exact solution of the adiabatic equations of 
motion Eqs. (|16p for real O c ff , as can be verified by 
insertion. If the integral appearing in the trigonometric 
functions is equal to it (a "7r-pulse" ) , all of the population 
ends up in the final state; if it is tt/2 (a "7r/2-pulse"), half 
of it. 

While operation in the Bragg regime is lossless, it re- 
quires relatively long interaction times. In the previous 
section, we used the condition 4(1 — n)u r 3> ft, which 
translates into 



l) n uj r 



2 n {n- l)! 2 



(19) 



This is ~ uj r for n < 5, but drops rapidly, e.g., f2 c ff "C 
2 x KT'V for n = 10. 

For later use, we consider the case of complex f2 c ff = 
|^off|e I¥ ', where <p is the argument of fi ff. Hermitic- 
ity then requires us to use f2* ff in the second of Eqs. 
(UHJ). The solution of these equations for constant f2 
is g- n (t) = cos 7j\Q e s\t, g n (t) = -ie~ lv sin ±\n eS \t. For 
time-varying complex fl, the generalization of Eqs. (fl8|) 
by substituting |f2 e ff| and inserting the factor of e~ lv into 
g n {t) is an exact solution for constant ip, and remains 
approximately valid as \(p\ <C |Q e ff|/2, as can be seen 
by inserting into Eqs. (fTB"]) . Thus, this solution holds 
for varying complex Q c s as long as its argument changes 
adiabatically. In the following, Q e g is always understood 
to be the absolute value unless otherwise stated. 



of the form 



to obtain 



g(z,t) = gt(t)g(z) 



(20) 



ig't ~ ut9t = o 
-^g" + 2{lcos 2 (kz)g-Lu t g = (21) 

with tot being the separation constant. The second equa- 
tion is the Mathieu equation. By using v = z/k, 



(jj r U)r 

n 

2w r ' 



(22) 



it can be brought to the standard form 

g"(v) + [a-2qcos(2v)]g(v) = 0. (23) 

This is an eigenvalue equation in a. The eigenfunctions 
can be expressed by Fourier series 

oo 

se 2n +i(v,q) = ^2 B 2r+i sin ( 2r + i) 11 > 

r=0 
oo 

Se 2 „ +2 (w,g) = B 2r+2 Sin ( 2r + 2 ) V > 



r=0 



ce 2n (v,q) = ^A^ n) cos2 



r=0 



ce 2n+ i(v,q) = Y, A 2r n + + i c° s ( 2r + 1 ) v - ( 24 ) 



r=0 



The eigenvalues associated with the v— even functions ce 
are denoted a^, those associated with odd functions se are 
denoted b r [5l| . Insertion of these series into the Math- 
ieu equation allows to determine the eigenvalues and the 
Fourier coefficients [SJ HI]. This is tedious, but tables 
[5l| and standard numerical routines allow to find nu- 
merical values easily. 

In this paragraph, we use units with u> r = 1. For an 
atom initially in a pure momentum state with p = —nhk, 
we express the wave function tp(t = 0) as a series of 
Mathieu functions: 



(25) 



F. Mathieu equation approach 

For constant r2 , we can apply the method of separa- 
tion of variables to Eq. Q . We are looking for a solution 



This is possible, because the functions ce and se form 
a complete orthogonal set [5l[. The coefficients of the 
expansion are thus given by the Fourier coefficients of 
the Mathieu functions. Once the they are known, we 
can write down the amplitude of finding the atom with 



a momentum nt at a later time: 



lf = ±- I*" e- in ' v ^^ce m (M)e ,(a "' +2 ' )i 

\ m 



,2iqt 



^A^A^e^ + B^B^e' 



bmt 



(26) 



(The e~ m f v has a minus in the exponent because this is 
a reverse Fourier transform.) While in general this is a 
very complicated function of t, for low values of q only 
si™ 1 ' 1 and A^f 1 will be large. If we neglect all others, 

\9n i ] 2 ^\{\A^\ i + \B^ 
±2\A^\ 2 \\B^\ 2 cos(a n ^b nt )t) . (27) 



The plus sign is for n/ = rij, the minus sign for n/ = — n,-. 
Thus, the atoms will oscillate between rij and — n, with 
an effective Rabi frequency f2 e ff = a ni — b ni . 

For an explicit example, let the two photon Rabi fre- 
quency be £1 = 2q = 3 for < t < T and otherwise. 
Suppose further that for t < the atom is in an initial 
state ip = e~ l3v = cos(3i> ) — i sin(3w) having momentum 
p = —3hk. The coefficients are most easily obtained by 
numerical calculation of the Fourier integral. 

The amplitude of finding the atom with a momentum 
+3 at a later time is given by Eq. 



(2n+l) |2 e ia 2 „+it 



\Bi 2n+1) \ 2 e lb2 



For the population in the initial state, we find 



(28) 



_|^( 2n + 1 )|2 e 46 2 „ + itj 



(29) 



1. Losses 

The solution (Fig. []} oscillates quickly around the 
mean as given by Eq. (|27|) . The frequencies of oscil- 
lation are relatively large compared to the effective Rabi 
frequency and depend on the pulse amplitudes. Thus, ob- 
serving them requires very accurate timing and control 
over the pulse amplitudes. Especially, because of inter- 
ference fringes in optical setups, it is hard to achieve an 
amplitude stability better than about 1%. Thus, these 
oscillations may be hard to observe in practice and the 
transfer efficiency very close to one at some of their peaks 
is not useful in practice. Most of the time, the population 
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FIG. 1: Population in the final state |<7n(t)| 2 as obtained from 
the Mathieu equation approach plotted versus time t in units 
of . 



that can practically be achieved will thus be more close 
to the mean value as given by Eq. (|27p . In our previous 
example, this reaches a maximum of 0.914, i.e., 8.6% of 
the population are not transferred to the final state. 



2. Phase shifts 

In this picture, the initial and final momentum states 
have the same energy. If they were freely propagating (no 
interactions with neighbor states), their wave function 
should thus exhibit the same phase factor exp(iEt/H). 
The interactions, however, cause a difference of the phase 
which can be seen by considering the ratio of the ampli- 
tudes 



tan (A = 



g-3ft) 



(30) 



which can be calculated in a straightforward way. The 
phase is most conveniently discussed in terms of <fi— tt/2, 
by subtracting the phase of 7r/2 which is expected in 
the pure Bragg regime, compare Eq. (| 1 T[) . As shown in 
Fig. [21 this is an oscillating function of time. Around 
the time for a 7r/2-pulse, the amplitude of the oscillation 
of the relative phase changes are given by tan (j> w <f> rs 
(A^) 2 m (B^) 2 . However, since the phase is an oscil- 
latory function, there are instances where the phase is 
larger or vanishes exactly. 

For the practical application of this in atom interfer- 
ometry, where this relative phase adds to the phase to be 
measured, two remarks of caution are appropriate: (i) As 
mentioned before, making use of the theoretically exact 
vanishing of the phase shift at particular times requires 
very accurate timing and control over the pulse ampli- 
tudes, (ii) In certain interferometer geometries, equal 
parasitic phases of subsequent beam splitters cancel out. 
However, since this depends sensitively on very small am- 
plitude changes of the pulses (that would affect the times 



7 




where Q is the vector g with g± n removed, Ti is TL with 
rows and columns of g± n removed, and 

c = • • 0, o 3 _„, n* 5 _ n , • • • o, n 5+n , sr g+n , • • • )• 

(34) 

Suppose for now that Q adiabatically follows g± n , i-e., 

iQ ps 0. C/ thus can be expressed as functions that are 
linear in c, and thus g± n (but not necessarily f2): 



FIG. 2: Phase cf> — n/2 of the final state relative to the initial 
state plotted versus time t in units of lo^ 1 ■ 



of the zero-crossings of the wiggles in Fig. [2]) , the cancel- 
lation is impaired. Thus, in practice, it may be impossi- 
ble to rely on this exact vanishing or cancellation. 



II. EFFICIENT METHOD FOR CALCULATING 
THE EFFECTIVE RABI FREQUENCY 

The adiabatic elimination process yields a simple equa- 
tion for the effective Rabi frequency Q e ff- However, it 
is inappropriate in the quasi-Bragg regime. Higher or- 
der corrections, that tend to reduce f2 e ff, will have to be 
taken into account. In this section, we shall determine 
these corrections for arbitrary scattering orders. A cal- 
culation for second order scattering has been published 
previously in [47| . 

The natural approach to determine Q e g is via the 
eigenvalues in the Mathieu equation formalism. This ap- 
proach can yield f2 e ff to any desired accuracy (for con- 
stant fi) by calculating the eigenvalues of the matrix rep- 
resenting an appropriately large subset of the infinite set 
of equations. In this section, we are looking for an effi- 
cient iterative method to calculate f2 G ff, for constant as 
well as time- varying f2. This is, at the same time, a 
method for calculating the difference of the eigenvalues 
a m — b m of the Mathieu equation. This method is based 
on an extension of the idea of adiabatic elimination. 

We express the equation of motion Eq (fT2|) as the ma- 
trix equation 



Ug = ig 

where the vector g contains the g n 

<?=("■ ,9-n,9-n+2, ■■■)■ 



(31) 



(32) 



In analogy to the above solution in the Mathieu equation 
formalism, we are interested in a solution that is slowly 
varying in time, in which the population is mainly con- 
sisting of <?_„ and g+n . The evolution of the other states 
is governed by the equation 



HQ = c + iQ, 



(33) 



Q H 1 ( iC, 

Let us define 



'(0) 



= ^(o) (ft, ft*) g _ n + £(o) ( ft ; ft*) g+ni 

that is, 
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(0) 

-n+2m 



d<$ (n,n*)g_ n + EW (n,n*)g+, 



(35) 
(36) 
(37) 



In 



ig-n 



(0) 
9-n-1 



o (0) 

y-n+2 



,(°) 



we apply Eq. ([57)) to replace the y_„ ±2 
analogy to Eq. (fTrJ)) , 



(38) 
and obtain the 



ig- n = fi ac (Q, Cl*)g- n + ~fi e ff (fy O*) g +n , (39) 



where 



(0) 



(0) 



(40) 

An analogous computation leads to the expansion for 

9+n, 



ig-r 



Q&cfJ—r 



a. 



-g+r, 



9. 



off . 



in analogy to Eqs. (fTB|) . where |f2 e ff| is the effective Rabi 
frequency. The leading order in Vt of 91 c q obtained this 
way is identical to the one obtained in previous section, 
see Eq. (|T5|) . 

Although initially the fast varying iQ was set to zero by 
the adiabaticity assumption, we now take into account a 
slowly varying part due to the adiabatic following, which 
has similar time scale as the initial and final states and 
thus cannot be ignored in Eq. (|35[) . In the remainder of 

this section, iQ refers to this slowly varying part only. To 

first order, iQ can be approximated from Q^: 



«- 4 |(flWj. ,41) 
Inserting into Eq. (|35p . 



(42) 
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Since is a function of fl and g± n , the time derivatives 
ig±n = ^rS±n-2 + %g±n+2 contain as well as still 
unknown corrections of Q. Therefore, care must be taken 
to properly separate various orders of corrections. We 

expand Q and Q as Q = + cf« + • • ■ and Q = + 
+ • • ■ , where ((5^) are an order of magnitude 
larger than Q^ l+1 ^ (<5( l+1 *>) and are functions of g± n and 



+EV\n,n,...) g+n . 



(43) 



We expand 



^ = E ( ^ ( Q *^- 2 + ^ +2 ) + *<7-«£ W 



l=V / -/n 



q=p-l 

+ ^ (fi* fl £>_ 2 + 0<^ +2 ) + ^g + J^ (44) 



Eq. (f3"3")) thus becomes 



p=0 



(p) 



p=0 



(45) 



Since the time derivative decreases one order of magni 
tude for each increase in p, the (p+ l)-th order in Eq. (l45l 
is 



HQ 



'(p+i) 



i(5 



(p) 



HQ 



'(0) 



(46) 



Thus, 



/+n- 



(47) 



We are now ready to describe an iterative procedure to 
obtain Q to any desired order: We start from as 
defined in Eq. (|3"5|) . Each component is known as a 
linear combination of g± n , and therefore D^ ',E^°> are 
also known. 

Now suppose that g<>i\DW,EW 

are known for q < p. 
In the last Eq. (|4"4"|) , each 9_n± 2 >3n± 2 can De expressed 
as a component of which, by Eq. (|4"3")> . is known 
as a linear combination of g± n . We insert this into Eq. 
(147)) to obtain again as a linear combination of 

<7± n . The coefficients of this linear combination are the 
#(p+i). The process can now be iterated to ob- 
tain the next higher order. 



The efficiency of this method is, in part, due to the fact 
that each order can be computed using the same inverse 
matrix Ti . (Moreover, fi. is a tridiagonal matrix, which 
helps in computing the inverse.) O e g to M-th order is 
thus obtained by plugging ^2 p=0 G^ into Eq.([35|) and 
finding the coefficient of g+ n as in Eq. (|39)) . 

Note that Ti and TL in principle are infinite- 
dimensional. However, for obtaining the effective Rabi 
frequency to order O (fl 2k ), it is sufficient to include the 
initial and final states, states in between, and k nearest 
neighboring states on each side, i.e., m = — k, ■ ■ ■ ,n + k. 
Including more states yields the same result. 

With this method, we explicitly calculate the effective 
Rabi frequency for Bragg diffraction orders of n < 29. 
They can be given as power series in Q and Q: 



a, 



UJ r 8 n_1 [(n — l)!] 2 \U; 



!-E a « 



(2j) 



3=1 



101 



2j 




(48) 



At first, this results in a list of numerical values for the 
coefficients a and j3 for each n. However, closed expres- 
sions as function of n can be found, which are listed in 
appendix [XJ 

Eq. (|4"5)) also allows us to give validity conditions for the 
simple adiabatic elimination method presented in Sec. 
HE 21 For this to be a good approximation, the correc- 
tions should be much less than 1. Thus, we obtain 



2 4 (n 2 



Hi «i, 



ci 

UJ r fl 



< 1. (49) 



For large n, the first of these conditions translates into 
fl <C 4n 3 / 2 w r , which is actually larger than the one given 
by the adiabaticity condition uj- n +2m > (n — l)uj r 51. 



Population in other states 

Summing up the population of states other than g± n 
at the end of a pulse (where g_„ ~ 0, g +n ~ 1), we obtain 



g-g* 



1 



2 5 (n 2 - l) 2 

(n 4 + 6n 2 + 1) 
+ 2 9 (n 2 -l) 4 



CI 



(50) 



The population lost into other states after the pulse is 
switched off, when fl = (l = 0, vanishes. The method 
presented in this chapter is not suitable for obtaining 
those losses, because the states other than g± n have been 
assumed to adiabatically follow their neighbors, and the 
losses are a non-adiabatic phenomenon. 
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However, as long as \g m ^±n\ 2 -C |.g±n| 2 , the effect of 
the losses on g± n and thus Q c g can be neglected. 

This iterative method, although powerful for calculat- 
ing f2 c ff, does not approach an exact solution. It can be 
seen from Eq. (j39|) that this method gives no wiggles in 
the sinusoidal change of the initial or final state popula- 
tion, while there are fast variations in the exact solution 
of a square pulse as shown in Mathieu function section. 
However, it approaches the solution for the initial and 
final state as averaged over the high-frequency wiggles, 
Eq. ([27]) and thus predicts the correct effective Rabi fre- 
quency. 



III. ADIABATIC EXPANSION 



A. Losses 

We are mainly interested in the population g n = 
g n (oo). For that, we can take out a phase factor and 
note that 9(t — t ) = 1. The absolute squares of these 
give the population in the neighboring states. These are 
closed, analytic expressions for the losses arising in the 
second order. As an example for a third order correction, 

x6(t 3 -t 2 ) cos (I f neff(ti)dtiV (54) 



To investigate the losses, corrections to the adiabatic 
method must be calculated. We will relabel the results 
g± n of the adiabatic method as They represent the 
first order adiabatic approximation. We now want to 
calculate corrections to the population of the states, 



7 (2) 



(51) 



where to first order only the initial and final state are 
nonzero. For calculating the second order, we insert g±n 

into Eqs. (fT2]) . Inserting from Eq. (fT7|) . we obtain 
population of the levels next to the initial and final states 
to second order 

ig-n±2 = 4(1 =F n)aj r g^ n±2 

+ iocos Q J il eS (l/)dt 



Wn±2 = 4(1 ± n)uj r g n±2 
2 



fisinQy n eS (t')dtj . (52) 



(n > 2). These are all states for which g' 2 ' ^ 0. The 

process can be iterated: From the g±„± 2 , corrections g^ 
can be obtained, and so forth. Here and throughout, 
we shall drop the superscript ( 2 ) as long as no confusion 
arises. 

These inhomogenous equations can be solved by stan- 
dard methods, such as variation of the constant or a 
Green's function: 



S-n±a(*) = ~J dt fl(t )cosQ J n e s(t')dt 



Xe -4i^n)"At-t o ) ( t _ t ^ 

0»±a(*) = ~\j dt Cl(t )sin(^ J° n cB (t')dt' 
Xe -4i(i± n ) Ur (t-t )^ t _ to y 



(53) 



B. Phase shifts 



(2) (3) 

To 3rd order, both g„± 2 contribute to g„ ■ 



7 (3) = ll 



t rt 3 

dt 3 n(t 3 ) / dt 2 n(t 2 ) 



r e -ii(l-n)u r (t 3 -t 2 ) , e -4i(l+n)w r (t3-t 2 )j 



Qeff(h)dh (55) 



X COS - 

2 

[The 6(t 3 ~t 2 ) has been omitted because it is one through- 
out the t 2 integration range.] Since this is a complex 
number, the phase of g^ + g^ will be shifted by some 
A<fi relative to the adiabatic result, g„ . 

An upper limit on A(£> can be derived as follows: We 
will assume that all the losses derived in the previous 
section in second order contribute to 100% to g\ . Fur- 
thermore, we will assume that the phase of these contri- 
butions is 90° shifted relative to g„ ■ This will result in 
an upper limit on Acf>, since in reality only a fraction of 
the losses in second order will be re-introduced into the 
final state to third order, and since a 90° phase of these 
contributions is the worst case. Then, 



|A0| < 



,(3) I 



l<?i 3) l<V7, 



(56) 



\gn + gh 

where the total losses are denoted £. 

The equations derived in this chapter are closed ex- 
pressions for the calculation of losses and phase-shifts for 
realistic, time-dependent envelope functions. We will ap- 
ply them to square and Gaussian pulses in the next two 
chapters. 



IV. SQUARE PULSES 

These expressions for the second order are easily solved 
for square pulses having a peak two-photon Rabi fre- 
quency of for < t < T and zero otherwise. For the 
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effective Rabi frequency, we can use the result from the 
simple adiabatic elimination method, because, as it will 
turn out, losses become very large already in the region 
where this is still valid. In this case, 



U n 2 '-j J dt cos ^— ) e 



\ ±i(l T n)u r ta 



9n±2 



2 Jo 



dto sin 



2T 



4i(l±n)a; r to 



(57) 



The integrals can be calculated in a straightforward way; 
the result can be simplified by considering the terms lin- 
ear in f2 only (noting that (w r T) _1 is proportional to 
Cl n /u)™)- This does not lead to reduced accuracy, since 
in this second order adiabatic expansion, terms of higher 
order can be assumed to be zero. We obtain 



9-n±2 = 



n i 

w r 8(1 =F n) 
. Q 1 



9n±2 = I 



ui r 8(1 ± n) 



3 4i(l±n)oi r T 



(58) 



The losses in the final and initial state are t = |g_ n _ 2 | 2 + 
|<?_„ +2 | 2 + | 5 „_ 2 | 2 + | 5n+2 | 2 = 2(|«? n _ 2 | 2 + | 5n+2 | 2 ). We 
obtain 



in 2 n 2 + l 
16^2 (n 2 -l) 2 ' 



(59) 



as plotted in Fig. [3] As an example, for n = 3, f2 = 3 
we have I — 0.088, in excellent agreement with the exact 
solution in terms of Mathieu functions (0.086). 

The square-root of this is at the same time an upper 
limit on the phase shift for a 7r/2 pulse, and indeed the 
exact solution as plotted in Fig. [5] verifies this. Square 
pulses are thus, unfortunately, not suitable for quasi- 
Bragg scattering, because f2 must be reduced strongly 
if low losses and phase shifts are desirable. Already the 
above example, with 9% losses and thus ~ 0.3 rad max- 
imum phase requires a 7r-pulse time of about 30/w r , or 
about 3 ms for cesium atoms. Reduction of this to below 
10~ 4 , as required for high-precision atom interferometers, 
would take a pulse time of about 10 6 seconds. 

In the next chapter, we will study pulses with a smooth 
envelope function, Gaussian pulses. These will exhibit 
much lower losses and phase-shifts for a given pulse du- 
ration. 



V. GAUSSIAN PULSES 

In this section, we specialize to Gaussian pulses: 

A = (le^ 2 '^ , (60) 
so that, from Eq. Q15p . 



^eff = ^™ I -z — 



i V 1 - 1 i 



8uj r ) [{n — 1) 



112 



(61) 




FIG. 3: Losses for square pulses for n = 2, ... 10 plotted 
versus Q/cu r . 



The case of Gaussian pulses is more difficult, (i) One the 
one hand, this is because the integrals are much harder, 
(ii) With Gaussian pulses, higher Rabi frequencies can 
be used. Thus, we need to take into account higher order 
corrections for the effective Rabi frequency. 

For now specializing to 7r-pulses, we have the condition 



- / n cS (t>)dt> = ii- 



i 

8u> 7 



V27R7 2 



[(n-1)!] 2 ^ 2' 
(62) 

that we will use to determine a. We insert into the ex- 
pressions for the amplitudes of the neighbor states and 
obtain 



9-n±2 



— oo 
7T 



9n±2 



x sin 



4 6rP V V 2^ 2 



(63) 



(64) 



— oo 
7T 



-erp 



A" c VV 2d 2 

We denoted erp(x) = l + $(x) and $ is the error function 



' dt. 



(65) 



For 7r/2 pulses, the factors of 7r/4 in the last three equa- 
tions will be replaced by tt/8. No simpler form of this 
integrals has been found. However, for 7r-pulses, we can 
use 



sin 



cos x + sin x 



V2 



V4 



(66) 



and $(a;) = — $(— cc) to obtain g- n ±2 = i<]n^2- 
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1. Inclusion of higher- order corrections to the effective 
Rabi frequency 

For Gaussian pulses, the peak effective Rabi frequency 
can reach a level beyond the region of validity of Eq. (fT5|) 
before the losses as predicted from the results of the last 
section become noticeable. In this region, however, the 
second order adiabatic solution is still a good approxima- 
tion, as long as the losses it predicts are still low (because 
then, the third order will be even lower). Therefore, we 
can extend the region of validity into the interesting re- 
gion, where the losses just start, by including higher order 
corrections to fl e g as calculated in section |TTJ In princi- 
ple, of course, adiabatic expansion to higher and higher 
orders is a systematic way of obtaining results of simi- 
lar and higher accuracy and simultaneously predict the 
changes in the effective Rabi frequency. However, this is 
very inconvenient method for Gaussian pulses, for which 
the integral is hard even in the second order. 

For a 7r pulse we have the condition 



1 f°° 7T 

- J n eS (t')dt' = - 



(67) 



We insert Eq. (|48|) and note that while the a coefficients 
and ^ are real, fin is imaginary. Cl c g is given by 
the absolute magnitude of the series and we expand the 
absolute value as a Taylor series in terms of the imaginary 
part that we truncate after the (Ct) 2 terms. 



L0 r 



B«-l[( n _l)!]2 ^ 
2" 



where 



^ = -M 1} | 2 +/3i 2) 



(68) 



(69) 



This approximation will be justified in retrospect, as we 
will find that transitions having low loss satisfy the con- 



dition Eq. (|74p . which gives a lower limit on a [53|. At 

this limit, the leading neglected terms (/3n and fior n ) 
are starting to become non-negligible and the compari- 
son to a numerical simulation shows an about 10% error 
in the predicted sigma, see Sec. IV Dl The inclusion of 
more terms would reduce the error here, but would lead 
to extremely lengthy expressions later. In the region of 
low loss, they are negligible. 

The condition for a 7r-pulse reads, 



1-/3- 



1 



-a 



(4) 



n 4 



v(6) 



1 

8uj r 



fi 6 



V2t 



[(n-1)!] 2 ^ 

o 2 



^7 y »• + 4 "" uJ b r \j n + Q "7 "2 
From this, a can be computed without difficulty, see Fig. 
ID In the following, we take into account the higher- 
order corrections up to order f2 n+6 and l/(uj r a) 2 . For 
the integrals appearing in the trigonometric functions, 
we find: 



2 Q ^ dt> = (8c r )«-M(n-l)!] 2 V 2n 




(71) 



Making the same assumptions as above, we can neglect 
the l/cr 3 term. Using Eq. ([70]) . we can write 



sin 



J 2 ) ^ rri ( U+2 t 



n+4'- 



1 

2 

, (4) n 4 



n eS (t')dt' 



1-/3—4- 



„(4) Q4 



The functions erp(yn/2£), erp(^/ (n + 2)/2t), . . ., al- though similar, show an increasing slope with respect 



12 




100 120 
fi[w r ] 



FIG. 4: a in units of uj t for 7r pulses for Bragg diffractions of 
order n = 4, 5, 6 (above) n = 7, 8, 9 (below) plotted versus fi. 
The effective Rabi frequency has been inserted up to order 

n n+6 . 



to t. As they are integrated over, and as y/(n + 2)/2 w 
yJn/2, this is insignificant for large n and we can take 
these functions to be equal. We then recover Eq. (|6"3")) . 
the only difference being the replacement of a by the so- 
lution of Eq. ([70|) . It remains to actually calculate the 
integral. 

This integral is a function of n and a. a, in turn, 
is determined by f2 for tt or tt/2 pulses. Unfortunately, 
the integral cannot be solved exactly. The general struc- 
ture of the integrand is a periodic oscillating function 
exp[— iu> r (l ± n)t] times an envelope that is peaked near 
t = 0. For g n ±2, the peak is limited on the right mainly 
by the Gaussian exp[— t 2 /(2a 2 )] and on the left by the 
sine of the error function. For low Rabi frequencies, a 
is high and thus the peak is broad compared to the pe- 
riod l/[(n± l)tu r ] of the oscillating function, which thus 
averages out to a very low value. When the peak width 
becomes comparable to the period, however, this is no 
longer the case and the value of the integral will increase. 
Since we are interested in the region where the losses are 
nonzero, but still low, it is sufficient to solve the integral 
with a method that is valid in this region. The method 
of steepest descent, or "saddle-point" method is suitable. 
As derived in Appendix [Cl neglecting a phase factor that 
is of no consequence unless one wants to proceed to higher 
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cr \co, 



FIG. 5: Numerical calculation of the integral 
f™ oo exp[-t 2 /(2cr 2 )-iuJr(l-n)t]erp{(w/4) sm[^/nj2{tj 'a)]}dt 
as plotted versus a. Parameter: n — 3, 4, ... 10. a is in units 
of inverse recoil frequencies. 



orders, 



9n±2 




l±n) ; 



aV 2 



n l/3p J ~ l 9~n T 2 ' 

(72) 

see Eq. ([CTO]) . Here, T m 1.64874 and D is a factor of 
the order one that is given in Eq. (|C12[) and plotted in 

Fig. nni 

Alternatively, the integral can be computed numeri- 
cally, see Fig. [5] Accurate numerical computation, how- 
ever, becomes difficult when the losses are low, because 
then the integrand will make very many oscillations. The 
saddle point method, however, works well especially in 
this region. Comparison of both shows good agreement 
for the range we are interested in, where the losses are 
lower than about 10%. 

As shown in Appendix [Cl the results can be easily 
adapted to 7r/2 pulses if T is replaced by — 1.5043 
and D by Z? w / 2 , which is also plotted in Fig. [TO] Thus, 
the losses for 7r and 7r/2 pulses are essentially the same. 



Losses 



The losses in the final and initial state are £ — 

\g-n-2\ 2 + \9-n+2\ 2 + |<M- 2 | 2 + |<M+ 2 | 2 = 2(| ff „_ 2 | 2 + 

|ffn+2| 2 )- Since the terms with the plus sign in the expo- 
nential can be neglected, we obtain a simple form 



A{VtDaye 



2 -16(l-n) 2 CT 2 ^/(n 1/3 r) 



(73) 



For losses below the 10 2 level, we have the condition 
(the factors outside the exponential are of order one) 



auj r > 



V2rin(2): 
4(n - 1) 



>V6 



0.38- 



i 1 / 6 
n — 1 



(74) 



Since the losses decrease rapidly with longer cr, very low 
losses can be reached: For example, if the theoretical 
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FIG. 6: Losses for Gaussian 7r pulses for Bragg diffractions 
of order n = 4 to 10 plotted versus S7. The effective Rabi 
frequency has been inserted up to order Q n+6 . 



losses should be lower than 10 10 (which is then clearly 
negligible compared to losses due to technical influences), 



,1/6 



<JL) r > 1.5- 



(75) 



We can now compute the losses and pulse durations as 
functions of fi, see Figs. [B]andU] 



B. Phase shifts 

As discussed, the parasitic phase shifts in radians are at 
most equal to the square-root of the losses, as computed 
in the previous section. Since the losses are such a steep 
function of f2, slightly reducing f2 below the level where 
the losses become appreciable will essentially reduce the 
phase shifts to negligible levels. For example, if Eq. ([75]) 
is satisfied, these theoretical shifts will be below 1CP 5 . 

For Gaussian pulses with an appropriate choice of O, 
the theoretical losses considered here are so low that they 
are negligible in practice. The practical losses will then 
be dominated by issues such as single-photon excitation, 
finite laser beam size, finite size and temperature of the 
atomic cloud, wavcfront distortions, and other things. 
These losses are based on entirely different mechanisms 
and will therefore not necessarily lead to a phase shift as 
the theoretical losses calculated in this paper. 



C. Truncated Gaussians 

Any experimental realization of the Gaussian must be 
truncated somewhere. Here, we consider the modifica- 
tions of the above considerations that arise if the Gaus- 
sian is truncated on the right side (only) at t = r a. 
The integrations in Eq. (|63|) will now run from — oo to 



t. It is convenient to write the final state after a trun- 
cated Gaussian pulse as g— n ±2,n±2 + <?- n ±2n±2J where 
g-n±2,n±2 is given by Eq. (53} and 



9-n±2 



4 ran 1 ^ J 



(76) 



x cos \ ^erp 



2<7 2 



9n±2 — r, 



(77) 



sin < — erp 



2cr 2 



a 4i(l±n)avt 



Since t > cr, we can use &(t) ps 1 — e~ t /(y/ni) and 
approximate the trigonometric functions to leading order 
of the argument. We obtain 



n l2ira 2 

9-n±2 — *77 



n 



^ e -(n+l)t 2 /(2tT 2 )+4i(l^n)w r t 



9n±2 = ™ / dte-* l [2 ° ') e «(i±»K 



(78) 



For t > t ^> cr, the Gaussian e t2 /( 2cr2 ) (and a fortiori 
e -(™+ 1 ) t /( 2cr )) i s a very steep function of t, which is 
non-negligible only in a small region above r. In fact, we 
can assume that the other functions in the integral are 
constant and have the values they take at t = r. The 
integrals then reduce to error functions, which in turn 
can be replaced by the asymptotic expression for large 
arguments. We obtain 



9-n±2 



Jn±2 



4 y/n{n + 1)t 2 

no- 2 



e -(n+l)T 2 /(2cT 2 ) e 4 l (l T n)"r7 



e -T7(2cT 2 ) e 4i(l±n)cv7 



(79) 



Truncation on the left interchanges the expressions for 

9n±2 and 9-n±2- 

For an example, consider r = cry' —2 hx(rf) is chosen 
such that the truncation happens at a fraction r\ of the 
peak amplitude. Obviously, the smaller r), the smaller 
the effect of truncation, which also decreases for large n, 
(as expected, as the Gaussian to the power of n has a 
narrower peak). Consider, for example, 77 = 1/100 (i.e., 
r fa 3cr) and n = 5. Then, \gL n±2 \ = 3.9 x lO^Cla. 
That means, the effects of truncation can be reduced to 
negligible magnitude. 

In order to keep phase-shifts low, it is desirable to 
choose r such that 4w r r = 27r, 47r, . . ., in which case the 
amplitudes for the neighbor states are real, i.e., there is 
no extra phase shift due to the truncation. On the other 
hand, purposely truncating the pulses at a variable r thus 
provides a means of testing the influence of these phase 
shifts in experiments. 
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0.5 1 1.5 2 

FIG. 7: Simulated populations of the initial and final state 
for n = 5, H = 38 plotted versus a in units of inverse recoil 
frequencies. 



D. Accuracy of the estimates and higher-order 
effects 

The losses are very steep functions of Cl. Therefore, 
the error in the predicted loss for a given Rabi frequency 
is much larger than the error in the inverse function, the 
Rabi frequency for a given loss. To confirm these re- 
sults, we numerically integrate the Schrodinger equation 
in the momentum representation. We include sufficiently 
many "outer" momentum states until the results are es- 
sentially unaffected by including more (usually, 5 are suf- 
ficient). The simulation program itself is checked for con- 
stant Rabi frequency against the exact solution in terms 
of Mathieu functions. 

The result for Gaussian pulses with n = 5 and f2 = 
38w r is shown in Fig. [71 For a < 0.1/oy, most of the 
population is still left in the initial state, |g-s| 2 ~ 1. 
For slightly larger a, the population gets driven out of 
the initial state. However, the transfer efficiency to the 
final state g§ is very low, as expected in the Raman-Nath 
regime. When the pulse gets longer, the system enters 
the quasi-Bragg regime and the losses become quite low. 
At <j = 0.52/ctv, the population of the initial state has a 
minimum of |<7-5| 2 ~ 0.0035 whereas the population of 
the final state |(7s| 2 ~ 0.89. The remaining 11% of the 
population are lost into other states. At a ~ 0.58/oj r , 
the final state has a first maximum of |(7s| 2 « 0.94; at 
this time, |g-s| 2 ~ 0.029, i.e., about 6% are lost into 
other orders. Thus, in this regime, the minimum of the 
initial and the maximum of the final state do not coincide. 
For even larger a, the losses become negligible and the 
system performs Pendellosung oscillations just as in the 
Bragg regime, although the adiabaticity criterion is still 
violated. 

For comparison, the second order adiabatic theory pre- 
dicts about 6% losses for Cl = 35, rather than 38. The 
bulk of the 10% discrepancy is due to the remaining error 



in the calculation of a from Eq. ((70)) . Another prediction 
of the theory that is confirmed by the simulation is the 
symmetry of the losses g n ±i = — <7* „ T 2- 

To third order, a part of the population that is lost to 
the neighbor states to second order returns to the initial 
and final state in the third order. This has two main 
effects: First, it reduces the predicted losses. Second, 
since the returning population will have a different time 
dependence for the initial and final state, it can account 
for the difference of the a that gives maximum population 
in the final state versus the one that gives a minimum in 
the initial state. We will not consider this. 



E. Practical example 

Suppose we want to achieve In = 16-photon Bragg 
diffraction with the lowest possible time compatible with 
losses below 1%. Eq. (1741) gives a minimum a > 
O.OSo;,? 1 w 6 /is. This is substantially faster than al- 
lowed by the adiabaticity limit Eq. (fT9)) . which would 
require fl c ff -C 0.007a-v ~ 27r x 14 Hz, i.e., transition 
times on the order of 0.1s. For applications in pre- 
cision atom interferometry, where the unwanted phase 
shift, that is estimated as the square-root of the losses, 
must be low, we might want to have theoretical losses as 
low as 10~ 10 (even if then additional losses due to tech- 
nical reasons will be much larger) . This can be obtained 
with a > 0.3W" 1 « 24 fxs, according to Eq. (f75|) . For this 
case, we can read off a required peak two-photon Rabi 
frequencies of fl ~ 80uj r from Fig. 2J 

For an explicit example, consider cesium atoms driven 
on the Q 2 S 1/2 ,F = 4, m F = -> 6 2 P 3/2 (D2) transition 
which has a wavelength of 852 nm, with a detuning of 
10 GHz. The relevant data of that transition is found in 
Ref. (54J. To determine the two-photon Rabi frequency, 
we neglect the hyperfine splitting of the excited state 
(which is on the order of a few 100 MHz) and sum over 
the excited states \F' = 5,mp = 1), \F' = 4, rap = 
1), and \F' = 3,mp = 1). From summing the matrix 
elements, we obtain a two-photon Rabi frequency of = 
§(/// sa t)±(r 2 /<S), or n/uj r = 0.20//(mW/cm 2 ), where 
/sat = l.lmW/cm 2 . 

If we have, for example, laser beams with a Gaus- 
sian waist of wo = 0.6 cm, the intensity at the center 
is I = 2P/(ttwq), where P is the total power. Inserting, 
we obtain Cl/uj r = 0.35P/mW. Thus, a power of about 
250 mW is necessary to reach D,/u> r ~ 88. 



VI. DISCUSSION, SUMMARY, AND OUTLOOK 

We have given an analytic theory of quasi-Bragg scat- 
tering. This is the range of short interaction times for 
which the usual adiabatic theory of Bragg scattering 
breaks down. Thus, population shows up in scattering 
orders other than those allowed by the Bragg condition. 
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However, we find that this population can still be ex- 
tremely low provided that the interaction is switched on 
and off with a smooth envelope function. 

For calculating these effects, we introduced a new 
method for solving the Schrodinger equation, adiabatic 
expansion. The Schrodinger equation in momentum 
space is a coupled system of differential equations de- 
scribing the population in the different momentum states. 
A first approximate solution is obtained in the usual way 
by adiabatically eliminating the states between the ini- 
tial and the final state. We then re-insert this first-order 
solution into the Schrodinger equation. This results in 
inhomogenous differential equations for the momentum 
states next to the initial and final state. Thus, we obtain 
a second-order solution. While this process can be iter- 
ated, the second order is sufficient for our purposes. Un- 
like previous approaches, this method allows us to treat 
arbitrary diffraction orders and envelope functions. We 
expect that this method is useful for obtaining higher ac- 
curacy results in all situations that are usually treated by 
adiabatic elimination. 

We also present an efficient method for calculating the 
effective Rabi frequency, which is related to the eigen- 
values of Mathieu functions. Closed expressions are ob- 
tained for the Rabi frequency up to eighth-order correc- 
tions for arbitrary Bragg diffraction order. 

We treat diffraction with a scatterer having square 
and Gaussian envelope functions as examples. Square 
envelopes lead to a high loss of the population into unde- 
sired momentum states (as experimentally observed, e.g., 
by [H). Moreover, due to couplings of the desired output 
state with the loss channels, the output state becomes 
phase-shifted as a function of the interaction strength and 
-duration. On the other hand, Gaussian envelopes can 
lead to extremely low losses, and hence phase shifts, even 
for pulse times that substantially violate the adiabaticity 
criterion. The effects of truncation of the Gaussian to a 
finite waveform can be made negligible by suitably cho- 
sen truncation. Comparison to a numerical integration 
of the Schrodinger equation verifies that the second order 
is sufficient for our purposes. 

While these results are important in all situations were 
Bragg diffraction is applied (like acousto-optic modula- 
tors) , we focus on the context of atom interferometry. For 
example, some present high-precision atom interferome- 
ters demand phase errors A(f> below 10 _5 rad [35| . For 
Gaussian pulses, the total losses I are given by Eq. ([75]) 
and A0 ~ y/£ < 1CT 5 can be satisfied by operating with 
a a longer than about 1.5/ui r , see Eq. (|75p . The strong 
dependence of i on f2 means that a slight reduction of fl 
corresponds to a strong reduction of A<p. On the other 
hand, for square pulses, \fl < 10 ~ 5 can only be fulfilled 
for f2 <C uj r , which will in practice mean unrealistically 
long transition times. 

While atom interferometry is a field to which our cal- 
culations are, we hope, useful, in this work we restricted 
attention to one single beam splitter. The application of 
the results to a full interferometer, e.g., in the Ramsey- 



Borde or other geometries, is a matter that we did not 
consider. While we considered the phase shift arising 
within one beam splitter, the signal of an atom interfer- 
ometer is given by the total amplitude of all interfering 
momentum states. In certain situations, the neighbor- 
ing momentum states (considered as "losses" in here) 
contribute to this, leading to an apparent distortion of 
the interference fringes. This has been considered for 
the Raman-Nath regime in [HBL l56| . An upper limit for 
this can be given by y/I, by assuming that all the lost 
population interferes and conspires to produce the maxi- 
mum effect. Thus, operation in the quasi-Bragg regime is 
suitable for reducing this contribution. Moreover, an ac- 
tual atom interferometer consists of more than one beam 
splitter. Since the undesired output states of the first one 
have different momentum, they will not be addressed by 
the subsequent (velocity-selective) beam splitters. This 
strongly reduces this effect. A further reduction is possi- 
ble by choosing the geometry of the interferometer such 
that the undesired output states do not interfere. 
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APPENDIX A: COEFFICIENTS OF A POWER 
SERIES EXPANSION OF THE MATHIEU 
EIGENVALUES 

We start from the numerical values obtained in Sec. HT1 
and make a polynomial ansatz D n = (n + ki)(n + fea) . . . 

for the denominator of the . The coefficients ki are 
then obtained by factoring the denominators into primes. 
If, for example, the prime 7 first appears at n = 7 — 1 
and n = 7 + 1, we choose one of the ki equal to — 1 and 
one to +1. This leads to a guess for the denominator. 
Subsequently, an ansatz is made for the numerators N n 
by a polynomial of Kth order. The coefficients are de- 
termined by solving X^fcLo a k nk = N n for the coefficients 
a/j. The resulting expression is then confirmed by direct 
comparison to the above calculation. As it turns out, 

(2?) 

the a>n , n = 1,2,.. . N can be given by polynomials of 
relatively low order compared to N. This suggests that 
the polynomials may turn out to be exact expressions for 
all n, although we have made no attempt to proof this 
conjecture. 

The functional forms of /3fp and fiffl are obtained as 
follows. For pn\ a factorial denominator (n — 1)! is 
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found, and then a recursive relation of the numerators is 
observed. After simplifying the relation, the functional 

(2) 

form is obtained and tested with larger n s. For p„ , 
u-n+m in fi are kept not evaluated, and a pattern of the 

(2) 

appearance of the u s in p„ is found. The pattern is then 
simplified with w's evaluated as w_ n +2 m = m(n — m)u) r 
thus yields the form. 

As a result, we obtain the following closed expressions 
for the coefficients: 



v(2) - 



„(4) 



-11141 



2 4 (n 2 - l) 2 ' 



7077888' 



«n>a = -(™ + 4)(4n 5 - I5n 4 
+64n+ lll)/[2 n (n 2 - 



-32n d 
1)V 



12n 2 
-4) 2 ]. 



(6) 



a 



1086647 
9555148800' 



(6) 



-872713 



1087163596800' 



°4>3 = ( n + 6)(4n 10 - 45n 9 + 76n 8 + 846n 7 
+484n 6 - 3960n 5 - 14824n 4 - 3078n 3 
+31904n 2 + 26973n + 30740) 
/[3! 2 14 (n 2 - l) 6 (n 2 - A) 2 (n 2 - 9) 2 ], 



APPENDIX B: MEAN OF THE EIGENVALUES 
a n ,b n OF MATHIEU FUNCTIONS 



f2 ac , that was introduced in Eq. (|39[) . is the common 
energy shift of the effective two-level system, which thus 
corresponds to the mean (a n + b n )/2 of characteristic 
values for real and constant £1. Thus, by taking into 
account of the difference in energy reference in Eq.Q 
and Ti., the mean of a n , b n can be obtained. As a result, 
both a n and b„ are known independently. 

In fact, for real constant O, every term in the power 
series expansion for the mean value as a function of n can 
be obtained efficiently The result is valid for all n as long 
as the denominator of the coefficients listed below do not 
vanish. This method is used to calculate the power se- 
ries expansion up to the 14th order in the parameter q as 
introduced in Eq. (|22[) on a personal computer in 10 min- 
utes, and the result agrees with Abramowitz and Stegun 
[5l[, who list the terms up to order q 6 . The subsequent 
terms are as follows: The coefficient of q 8 is 

(274748 + 827565n 2 + 64228n 4 - 140354n 6 + 9144n 8 
+1469n 10 )/[2 2 V 2 - l) 7 (n 2 - 4) 3 (n 2 - 9){n 2 - 16)]; 



of q 1 



,.(») 



..(«) 



,.(») 



-20778032863 

2254342434324480' 
16738435813 

2727476031651840000' 
-218963004049 

2301307901706240000000' 



a nL = ~(n + 8)(5191881936 + 5562082816™ 



l n>4 

+9349559664n 2 4 
-3794628570n 5 - 
+888115497n 8 - 
-109 6 6 5671 1 



953069376n 3 - 
375347622n 6 4 
■ 172767750?i 9 - 
1 + 16309509^ 



6361852029n 4 
1490778180?! 7 
■ 198150468n 10 
+ 3762906n 13 



-472854n 14 - 257532n 15 + 11487n 16 + 4998n 17 

-720?i 18 + 32n 19 ) 
/[4!2 21 (n 2 - l) 8 (n 2 - 4) 4 (n 2 - 9) 2 (n 2 - 16) 2 ], (Al) 



ra-1 



k=l 



= — 



71-2 



k-1 



k=2 



kin — k) 



# fe -i(A2) 



At present, we have proved these expressions for n < 30 
by direct computation, as explained above. The ex- 
pressions have also been checked against the differences 
o>n — b n of the eigenvalues of Mathieu equation, as far as 
they are available from [511 ]. 



(4453452 + 20651309n 2 + 13541915n 4 - 2844430™ 6 
-1039598n 8 4- 69361n 10 + 4471n 12 ) 
/[2 24 (n 2 - l) 9 (n 2 - 4) 3 (n 2 - 9)(n 2 - 16)(n 2 - 25)], 



of q 1 



(1155192131376+ 6474038960008n 2 
+4470328527807n 4 - 3667584923421n 6 
-518221243968n 8 + 604333473552n 10 
-86398056330n 12 - 5960020482n 14 
+2031809256n 16 - 119406048n 18 - 131341n 20 
+121191n 22 )/[2 30 (7i 2 - l) n (n 2 - 4) 5 (n 2 - 9) 3 
x(n 2 -16)(n 2 -25)(n 2 -36)], 

and of g 14 , 

(41218724372688 4- 302286294466120n 2 
+433922366490105n 4 - 32881683767026n 6 
-136265027585703n 8 4- 21249098173752n 10 
4-79681 13847666n 12 - 1857745713708n 14 
4-30568658442n 16 4- 15106984464n 18 
-880061323n 20 4- 3059598n 22 4- 441325™ 24 ) 
/[2 33 (n 2 -l) 13 (n 2 -4) 5 (n 2 -9) 3 
x (n 2 - I6)(n 2 - 25){n 2 - 36){n 2 - 49)]. 

This is the first time the close form higher order power 
series expansion terms of the mean characteristic values 
are reported. 
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APPENDIX C: CALCULATION OF THE 
INTEGRAL 

In this appendix, we calculate the integral 

rn t 



die- * /< 2CT ) sh 



4 erP r V 2a 



We substitute u^Jn/2{t/a), to bring the integral into the 
form 



9n±2 



12a 2 Q 

— g± 

n Iv 



where 



9± 



e -at 2 ^bt gm 



in Qerp(t)) 



dt 



(CI) 



(C2) 



and a — 1/n, b — —4(1 ± n)^/2a 2 /nuj r /v. A useful 
approximation is the saddle point method. We write 



9± 



f = —at 2 + In ^sin — erp(t) 



(C3) 



We expand into a Taylor Series near the maximum of / 
at to i 



wi 1 
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FIG. 8: to for 7r (lower graph) as well as n/2 (upper graph) 
pulses, plotted versus 1/a. 
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1.55 
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C 



.'/ : / '^{f(to) + \f"{t )(t-t ) 2 + 



Thus, 



9 



2vr o/(to) _ l&to+ _^_ 



ibtj dt . 
(C4) 

(C5) 



, nto) 

For finding to, we calculate 

/'(to) = -2at + ^e-'° cot [ Jerp(t )j - . (C6) 

to can be calculated numerically and is plotted versus 1 /a 
in Fig. E 

For the second derivative, wc calculate 



f"(t) = -2a- e"* V?rtcot -erp(t ) 



-21- 



csc 



rerp(to) 



(C7) 



To obtain a simplified expression for /"(to), we derive 
fromEq. (fC6|) : 



e cot 



esc 



erp(to) 
erp(t ) 



4at n 



/ 7T 



1 + ^a 2 t 2 e 2t ° . (C8) 

7T 



10 



15 20 
l/o 



FIG. 9: C (upper graph) and CV/2 (lower graph) plotted ver- 
sus 1/a. 



For obtaining the latter result, esc 2 x = 1 + cot 2 £ has 
been used. We obtain /"(to) = — 2Ca 2 / 3 where 



C = 



(a + 2a 2 t 2 + Je" 2t o + 2at 2 ) /a 2 / 3 



(C9) 



to as well as A depend on a only. As can be seen from 
Fig. [3J C is of order unity and does not vary strongly. 
Indeed, C has a maximum of T « 1.64874 at 1/a = 8.984 
and can be replaced by T with less than 3% error for 
3 < 1/a < 30. 

Inserting into the expression for the integral in the sad- 
dle point method yields 



.flD 2a 



7n±2 = -i \ exp<4i 

v V n 

2^.2, ,2 



(l±n) 2a 



U) r t(} 

n 



{l±n) 2 a 2 LJ 2 



~9n±2- 



(C10) 



va 2 /^nC 

The overall scaling is predominantly determined by the 
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The calculation of the previous sections can be carried 
through in analogy, inserting factors of 1/2. For example, 
equation (|63|) will have factors of 7r/8 in the cosine and 
sine functions, rather than 7r/4. This means, that g— n ±2 
will no longer be equal to —g* L ^2i so we have to treat 
these cases separately. In the following, we specialize on 
5n±2! the other ones can be treated in analogy. For the 
evaluation of the integral in the saddle point method, the 
changes can be summarized by replacing C and thus D 
by 
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FIG. 10: The factor D (upper graph) and D n / 2 (lower graph) 
plotted versus 1/a. 



factor 





b 2 - 




exp 


4Ca 2 / 3 _ 


= exp 



{l±n) 2 a 2 m 2 
va 2 / 3 nC 



(Cll) 



which is a strong function of a and, thus, Cl/ui r . The 
factor 



■ sin 



7T . . 



2 V Ca 2 / 3 
is plotted versus 1/a in Fig. [T0l 



e -at 



(C12) 



a. ty/2 pulses 
7r/2 pulses can be treated in analogy by requiring 

i f n cS (t)dt = J . (C13) 



C n/2 = {a + 2a 2 t 2 + —e- 2t °+2atl)/a 2 l\ 



D 



7T/2 



2yJC n/2 nW 



rerp(to) 



-tl/n 



(C14) 



(note that also to will have another value). CV/2 has a 
maximum value of ~ 1.5043 at n — 9.534 and can be 
replaced by that value with little error for 3 < n < 30. 

is also plotted in Fig. [TO] It is a slowly varying 
function of n with a value somewhat lower than that for 
7r pulses. However, in general losses for ir/2 pulses will be 
higher, because of the lower a that enters the exponential 
which chiefly determines the magnitude. 
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